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Abstract 

We report the first measurement of the forward-backward asymmetry and the ratios of Wilson 
coefficients A9/A7 and A10/A7 in B — ► K*£ + l~, where I represents an electron or a muon. We 
observe a large integrated forward-backward asymmetry with a significance of 3Aa. The results 
are obtained from a data sample containing 386 xlO 6 BB pairs that were collected on the T(45) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 11.30.Er, 11.30.Hv, 12.15.Ji, 13.20.He 
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Flavor-changing neutral current b — > s processes proceed via loop diagrams in the Stan- 
dard Model (SM). If additional diagrams with non-SM particles contribute to such processes, 
the decay rate and other properties are modified. Such contributions may change the Wil- 
son coefficients £0] that parametrize the strength of the short distance interactions. The 
b — > s£ + £~ amplitude is described by the effective Wilson coefficients Cj S , Cg S and Cf^ , 
whose terms have been calculated up to next-to-next-to-leading order (NNLO) in the 
SM. 

The magnitude of C| ff is strongly constrained from measurements of B — > X s ^ 0, 0] 
and a large area of the (Cf , Cf^ ) plane is excluded by branching fraction measurements 
of B -> K(*H+£- and B -> X s £ + £~ @, 0, 0, 0. However the sign of Cf and values of 
Cg S and C^q are not yet determined. Measurement of the forward-backward asymmetry 
and differential decay rate as functions of q 2 and 9 for B — > K*£ + £~ constrains the relative 
signs and magnitudes of these coefficients 0, Here, q 2 is the squared invariant mass 
of the dilepton system, and 9 is the angle between the momenta of the negative (positive) 
lepton and the B (B) meson in the dilepton rest frame. The forward-backward asymmetry 
is defined using the differential decay width, g(q 2 ,9) = d 2 T/dq 2 d cos 9 as 



2i j\ sgn(cos#)g(g 2 , 9) dcos9 



Li9{q 2 ,d)dcos9 

The numerator in Eq. ^does not cancel due to interference between the electroweak penguin 
and box diagrams, and can be expressed in terms of Wilson coefficients as 

l 

sgn(cos 9)g(q 2 , 9)d cos 9 

i 

= "CTo £(<? 2 ) {^<Cf)F x + ^CfF 2 ) , (2) 

where £ is a function of q 2 , and F lj2 are functions of form factors (the full expression can be 
found in Ref. HJ). 

In this Letter, we report the first measurement of the forward-backward asymmetry and 
ratios of Wilson coefficients in B — > K*£ + £~ . We use a 357 fb _1 data sample containing 386 
xlO 6 BB pairs taken at the T(AS) resonance. We also study the B + — > K + £ + £~ mode, 
which is expected to have very small forward-backward asymmetry even in the existence of 
new physics [12]. Charge-conjugate modes are included throughout this Letter. 

The data were taken at the KEKB collider [l3| and collected with the Belle detector [lij . 
The detector consists of a silicon vertex detector, a central drift chamber, aerogel Cherenkov 
counters, time-of-ffight scintillation counters, an electromagnetic calorimeter, and a muon 
identification system. 

The event reconstruction procedure is the same as described in our previous report 0. 
The following final states are used to reconstruct B candidates: K*°£ + £~, K* + £ + £~, and 
K + £ + £~, with subdecays K*° -> K+ir', K* + -> K° s ti + and K + tt°, K° s -> tt+vt", and 
7T° — > 77. Hereafter, K*°£ + £~ and K* + £ + £~ are combined and called K*£ + £~. 

We use two variables defined in the center-of-mass (CM) frame to select B candidates: 



the beam-energy constrained mass M bc = \J (E£ e&m / c 2 ) 2 — (p* B /c) 2 and the energy difference 
AE = E* B — E^ eam , where p* B and E* B are the measured CM momentum and energy of the 
B candidate, and E^ e&m is the CM beam energy. When multiple candidates are found in an 
event, we select the candidate with the smallest value of \AE\. 
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The dominant background consists of BB events where both B mesons decay semilep- 
tonically. We suppress this using missing energy and cos8 B , where 9* B is the angle between 
the flight direction of the B meson and the beam axis in the CM frame. These quantities 
are combined to form signal and background likelihoods, £ sig and C B bi an d event selection 
is then performed using the ratio 7Z B b = C sig /(C sig + C B b)- The continuum (e + e~ — > qq, 
q = u,d, s, c) background is suppressed using a likelihood ratio lZ rnri t (defined similarly 
to TZbb) that depends on three variables; a Fisher discriminant [15| calculated from the 
energy flow in 9 cones along the B candidate sphericity axis and the normalized second 
Fox- Wolfram moment [16j, the angle between the beam axis and the CM sphericity axis 
calculated with tracks used in the B meson reconstruction, and cosd B . Backgrounds from 
B — > J/ipX s , ip{2S)X s decays, below referred to as B —> ipX s , are rejected using the dilepton 
invariant mass. Backgrounds from photon conversions and 7r° Dalitz decays are suppressed 
by requiring the e + e~ invariant mass to be above 140 MeV/c 2 . 

The signal box is defined as \M^ C — < 8 MeV/c 2 for both lepton modes and 
—55 (—35) MeV < AE < 35MeV for the electron (muon) mode. We optimize the selections 
on lZ corit and 1Z B b f° r eacn K* decay mode and each lepton mode to maximize sensitivity 
to events with q 2 < 6 GeV 2 /c 2 . 

To determine the signal yield, we perform an unbinned maximum-likelihood fit to the 
Mbc distribution for events that lie within the AE signal window. The fit function includes 
signal, cross-feeds and other background components. The cross-feeds are misreconstructed 
}{(*)£+£- events with correct ("CF") and incorrect ("IF") flavor tagging. The cross-feed from 
X s £+£- events other than K^*H + t is negligible. The other backgrounds come from dilepton 
background, combinatorial K^'i^h^, K^h + h~ and ipX s events, where h represents a pion 
or a kaon. The dilepton background refers to the sum of all background sources with two 
leptons where the lepton is from (semi)leptonic meson decays, photon conversions and 7r° 
Dalitz decays. The K^h + h~ is from both combinatorial background and B meson decays. 

The shape for cross-feed events is parametrized by a sum of an ARGUS function |l7[ 
and a Gaussian whose parameters are determined from Monte Carlo (MC) samples. The 
dilepton background is characterized by an ARGUS function. The shape of each background 
is determined from a MC sample. (The K^e ± ji T background shape is found to be consistent 
in MC and data.) Since the shape for K^i^hr^ is similar to that for the dilepton background, 
we use the same parameterizations for both backgrounds. The residual background from 
ipX s is estimated from a MC sample of ip inclusive events and parametrized by the sum 
of an ARGUS function and a Gaussian. The background from events with misidentified 
leptons is also parametrized by the sum of an ARGUS function and a Gaussian. In the fit, 
all background fractions except the dilepton background are fixed while the signal fraction 
is allowed to float. 

Figure ^ shows the fit result. We obtain 113.6 ± 13.0 and 96.0 ± 12.0 signal events for 
K*£ + t and K + £ + e~, respectively. 

We use B — > K*£ + £~ candidates in the signal box to measure the normalized double 
differential decay width. For the evaluation of the Wilson coefficients, the NNLO Wilson 
coefficients C{ of Ref. |2j are used. Since the full NNLO calculation only exists for q 2 /m 2 < 
0.25 region, we adopt the so-called partial NNLO calculation |7[ for q 2 /m 2 > 0.25. The 
higher order terms in the C{ are fixed to the SM values while the leading terms A{, with 
the exception of A?, are allowed to float. Since the branching fraction measurement of 
B — > Xg'-f is consistent with the prediction within the SM, Aj is fixed at the SM value, 
—0.330, or the sign-flipped value, +0.330. We choose A 9 /A 7 and A w /A 7 as fit parameters. 
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FIG. 1: M bc distributions for (a) B -» K*£+£~ and (b) B -» samples. The solid and 

dashed curves are the fit results for the total and background contributions. 

The SM values for A 9 and A w are 4.069 and -4.213, respectively 0]. To extract the these 
ratios, we perform an unbinned maximum likelihood fit to the events in the signal box with 
a probability density function (PDF) that includes the normalized double differential decay 
width. The PDF used for the fit consists of terms describing the signal, cross-feeds and 
backgrounds: 

P(M bc , q 2 , cos 0; A 9 /A 7 , A 10 /A 7 ) 
= j-r- j si g e S i g (<? 2 , cos 0) g(q 2 , cos 0) 

1 "sig 

+ ^^icFecF(<? 2 ,cos%(g 2 ,cos6») 
iv CF 

+ -rr- fweiF (q 2 , cos 8) g(q 2 , -cos 9) 

iViF 

+ (1 — /sig — fcF — flF — fK'hh — fipXs) X 

{(fK*thPKHh(q 2 , COS 9) + (1 - fK*ih)V&\(q 2 , COS 0)} 

+ fK*hhV K *hh(q 2 , cos 9) + f 1 p Xs V^x s {q 2 ,cos9). (3) 

Here, VK*eh, Vdi, VK*hh and P^,x s are the probability density functions for K*£h, dilepton 
background, K*hh and ^X,; respectively. The quantities e S i g (iV sig ), ccf (^cf) and ejp 
(Nip) correspond to the efficiency function (normalization) of each signal and cross- feed 
component. Each fraction / is the probability of finding the corresponding component in the 
data sample for a given M bc value determined from the M bc fit, with the exception of f'K*m, 
which is the fraction within the dilepton background component determined from the MC 
samples.) The functions e and V for the dilepton background, K*£ and ipX s are obtained 
from MC samples. The K*h + h~ background shape VR*hh is obtained from K*h + h~ events 
and the momentum- and angular-dependent hadron to lepton misidentification probability. 

The renormalization scale /i is set to 2.5 GeV as suggested by Ref. 0|. The double 
differential decay width includes the form factor parameters and the bottom quark mass m&. 
We choose the form factor model of Ali et al. @,[T3| and a bottom quark mass of 4.8 GeV/c 2 . 

First, we measure the integrated asymmetry ^4fb, which is defined as 

j = I I-i sgn(cos 9)g(q 2 , 9)d cos 9dq 2 
FB IS-i9(q 2 ,8)d cos 9dq* ' 1 ' 
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We determine the yield in each q 2 and forward-backward regions from a fit to the M^ c 
distribution. Then we correct the efficiency and obtain 

A FB (B -> K*£ + r) = 0.50 ±0.15 ±0.02, 
A FB (B + -> K+£+r) = 0.10 ±0.14 ±0.01, (5) 

where the first error is statistical and the second is systematic. A large integrated asymmetry 
is observed for K*£ + £~ with a significance of 3.4c. The result for K + £ + £~ is consistent with 
zero as expected. 

We fit the K*£ + £~ candidates with the PDF of Eq. El The fit results of ratios of Wilson 
coefficients are summarized in Table Figure 121 shows the fit results projected onto the 
background-subtracted forward-backward asymmetry distribution in bins of q 2 . 

TABLE I: Ag/Aj and A\q/Aj fit results for negative and positive Aj values. The first error is 
statistical and the second is systematic. 





Negative A 7 


Positive Aj 


A 9 /A 7 
A 10 /A 7 


-15.3 ±H± 1.1 
10.3 ±|J ± 1.8 


-16.3 ±J;£ ±1.4 
11.1 ±H± 2.4 




FIG. 2: Fit result for the negative Aj solution (solid) projected onto the background subtracted 
forward-backward asymmetry, and forward-backward asymmetry curves for several input parame- 
ters, including the effects of efficiency; A7 positive case (Ay = 0.330, Ag = 4.069, Aiq = —4.213) 
(dashed), A\q positive case (A7 = —0.280, Ag = 2.419, A\q = 1.317) (dot-dashed) and both A7 and 
j4io positive case (A7 = 0.280, Ag = 2.219, Aiq = 3.817) (dotted) 9]. The new physics scenarios 
shown by the dot-dashed and dotted curves are excluded. 

We estimate contributions to the systematic error due to uncertainties in the physics 
parameters, finite q 2 resolution, efficiency and signal probability. We vary the value 
within the range allowed by the branching fraction of B — > X s j 0]. The bottom quark 
mass rrif, is varied by ±0.2 GeV/c 2 . The systematic uncertainty associated with the choice 
of the form factor model is taken from the difference in fit results using the models of Ali 
et al. and Melikhov et al. |l9j ]. The effect of q 2 resolution is estimated using a toy MC 
study. The effect due to cos# resolution is found to be negligible. The uncertainty in the 
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efficiency is estimated by changing the efficiency for pions with p < 0.3 GeV/c, electrons 
with p < 0.7 GeV/c and muons with p < 1 GeV/c by 10%, 5% and 10%, respectively, to 
obtain revised efficiency functions for signal and background PDFs. We change the shape 
parameters for the signal or background probability functions / and take the difference as an 
uncertainty in the signal fraction. The parameters are modified by ±ler for signal, dilepton 
background and K*h + h~ . We vary the normalization for cross-feed events and ipX s by 100% 
since we cannot determine the uncertainty from data. To assign the uncertainty in K*£ ± h T , 
we change the fraction fx* eh by 20%, which corresponds to the difference between MC and 
sideband events. Table HP summarizes the contributions to the systematic error. 



TABLE II: Summary of systematic errors. 



Source 


Negative A 7 


Positive A 7 




Ag/A 7 


Ai /A 7 


A 9 /A 7 


A 10 /A 7 


A 7 [18J 


+0.2 
-0.0 


±0.0 


+0.1 
-0.2 


+0.3 
-0.1 


m b (4.8 ± 0.2 GeV) 


±0.7 


±0.5 


±0.6 


±0.4 


Model dependence 


±0.7 


±1.7 


±1.0 


±2.2 


q 2 resolution 


±0.3 


±0.4 


±0.3 


±0.4 


Efficiency 


±0.1 


±0.0 


±0.1 


±0.1 


Signal probability 


+0.4 
-0.5 


+0.2 
-0.3 


+0.4 
-0.5 


±0.4 


Total 


±1.1 


±1.8 


+1.3 
-1.4 


+2.4 
-2.3 



The fit results are consistent with the SM values Ag/A 7 = —12.3 and Aiq/A 7 = 12.8. 
In Fig. El we show confidence level (CL) contours in the (Ag/A 7 , A 10 /A 7 ) plane based on 
the fit likelihood smeared by the systematic error, which is assumed to have a Gaussian 
distribution. We also calculate an interval in A 9 A 10 /A 7 at the 95% CL for the allowed A 7 
region, 

- 14.0 x 10 2 < A 9 A 10 /A 2 7 < -26.4. (6) 

From this, the sign of Agv4 10 must be negative, and the solutions in quadrants I and III of 
Fig.|3]are excluded at 98.2% confidence level. Since solutions in both quadrants II and IV are 
allowed, we cannot determine the sign of A 7 Aio. Figure El shows the comparison between the 
fit results for the negative A 7 value projected onto the forward-backward asymmetry, and 
the forward-backward asymmetry distributions for several input parameters. We exclude 
the new physics scenarios shown by the dotted and dot-dashed curves, which have a positive 
value. 

In summary, we have measured the ratios of Wilson coefficients in B — > K*£ + £~ decay 
for the first time by studying the forward-backward asymmetry in the angular distribution 
of leptons. We observe a large integrated forward-backward asymmetry with a significance 
of 3.4a. The fit results are consistent with the SM prediction and also with the case where 
the sign of A 7 A 10 is flipped. We exclude new physics scenarios with positive A 9 A 10 at 98.2% 
confidence. 
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We thank the KEKB group for excellent operation of the accelerator, the KEK cryogenics 
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FIG. 3: CL contours for negative A^. Curves show la to 5 a contours. The symbols show the fit 
(circle), SM (triangle), and ^lo-positive (star) cases. 
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